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Non-Perturbative Improvement of the Anisotropic Wilson QCD Action* 
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We describe the first steps in the extension of the Symanzik O(a) improvement program for Wilson- type quark 
qq ■ actions to anisotropic lattices, with a temporal lattice spacing smaller than the spatial one. This provides a fully 
relativistic and computationally efficient framework for the study of heavy quarks. We illustrate our method with 
^\ ■ accurate results for the quenched charmonium spectrum. 
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The simulation of charm or bottom quarks on 
(a sequence of) lattices with am q <g; 1 is currently 
not affordable, even in the quenched approxima- 
tion. Since standard lattice actions break down 
when am q > 1, one must design special actions to 
study heavy quarks. Currently two approaches 
are widely used: (i) Non-relativistic lattice QCD 
(NRQCD) Q, and (ii) the "heavy relativistic" or 
"Fermilab" approach 

The heavy relativistic approach aims to imple- 
ment the improvement program for the Wilson 
action at any quark mass. Then the action only 
has manifest 0(3) symmetry, and to improve it to 
a given order in the spatial momenta requires all 
coefficients to be mass-dependent. The approach 
we would like to advocate J|,|| can be denoted 
as "anisotropic relativistic", for short, ft is simi- 
lar to the heavy relativistic one in that the same 
terms have to be included in the lattice action, 
but it avoids the mass dependence of the coeffi- 
cients (except for one, see below) by using a tem- 
poral lattice spacing ag so small that a^m q <C 1, 
and aiming to improve the action fully only up to 
0(a). ft has become clear in the last few years 
that a non-perturbative determination of the lead- 
ing improvement coefficients in a Wilson-type ac- 
tion is very important . In our approach this 
can be achieved, we believe, but in others the 
large number and/or mass-dependence of the co- 
efficients make this appear very difficult. 

Consider an anisotropic lattice with fixed 
(renormalized) anisotropy £ = a/ao- Up to 0(a) 
a Wilson-type quark action on such a lattice can 



have the following terms, in continuum notation, 

f (l + O(am)), p(l + 0(am)), 

aD 2 , aD 2 ., ofp,^ ], aa kF ok , aa k iF kl (1) 

Here a^ v = -|[7 M ,7r/] and F^ v is the field 
strength. Using a field transformation we can set 
the coefficient of ft (say) to 1 and adjust the 
coefficients of the first three terms in the second 
line of (|j). We set the third term to zero, and ad- 
just the second-order derivatives to combine with 
the first-order ones to give the computationally 
efficient Wilson operators, = V M — ia M A M , 
where V M and are the standard discretizations 
of first- and second-order covariant derivatives 
(other possible strategies are discussed in 
This gives a quark matrix Q of the form 



Q = ni + Wo7o + v ^2 W k -yk 



c Q k F ok + lo 2^ v ki Fki 



k<l 



(2) 
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Here v is a "bare velocity of light" and u>o, oj are 
the temporal and spatial "clover coefficients" ; all 
three have to be tuned on the quantum level. 

It is instructive to compare our and the heavy 
relativistic approach. The latter corresponds to 
the special case of an isotropic lattice, £ = 1, and 
one must adjust the coefficients u,ujq,u} to all or- 
ders in the mass. As an example consider the 
classical tuning of v. Defining an "effective ve- 
locity of light" c(p) via the dispersion relation 
E(p) 2 = E(0) 2 + c(p) 2 p 2 , we demand c(0) = 1 
for a free quark. The result for v is shown in fig- 
ure [l] for various £, as a function of the pole mass 
m = ln(l + aomo)/ao of the quark. Note that 
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Figure 1. Classical bare velocity of light. 
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Figure 2. Dispersion relation -E(p) of a free, mas- 
sive quark for p oc (1, 1, 0) after tuning v. 

v has a strong mass dependence in the isotropic 
case, but only a very weak one for anisotropic 
lattices. The dispersion relation is shown in fig- 
ure ||. It is much better for anisotropic lattices, 
illustrating that they are the natural habitat for 
heavy quarks. 

One can also work out the exact classical mass 
dependence of loq and u). It is (almost) iden- 
tical to figure [I], with a strong/weak mass de- 
pendence in the isotropic/anisotropic cases. In 
our approach we can choose ujq and u> to be 
mass- independent, but must adjust v to 0(a) 
(for details and subtleties see 0). We can also 
tune v non-perturbatively, by demanding that the 
pseudo-scalar meson, say, satisfy c(0) = 1. In fact, 
it is simpler to determine v exactly for the heavy 
masses of interest in this manner; we thereby 
avoid expensive simulations at small masses to 
determine the 0(a°) and 0(a) pieces of v(m,£). 

The complete anisotropic Wilson action con- 



sists of the quark action described above and 
the anisotropic Wilson gauge action (see e.g. @), 
which is parameterized by (3 and the bare 
anisotropy £o- For given /3, £ and mo, tuning 
the full action involves determining £o> v, ujq and 
u> self-consistently to give the same renormalized 
anisotropy £ in the gauge and quark sectors, and 
no 0(a°,a) errors. Although not a problem of 
principle, in full QCD this is a costly tuning prob- 
lem, so we will here only consider the quenched 
case. Then the relation £ = £o can be de- 
termined once and for all in the pure gauge the- 
ory, which we did to high precision in Q| . 

To non-perturbatively tune the quark action 
one can use the following strategy: (i) Deter- 
mine v{m, £) by demanding that the pseudo- 
scalar have a relativistic dispersion relation for 
p — > 0. (ii) Determine the clover coefficients using 
the Schrodinger functional with the PCAC rela- 
tion as improvement condition (for details cf. 
|(|). For wo impose fixed boundary conditions 
in the time direction to generate a color-electric 
background field; similarly, for lu impose fixed 
boundary conditions in a spatial direction. 

On the classical level this strategy completely 
decouples the tuning of the three coefficients. On 
the quantum level this will not hold exactly, and 
one will have to tune the coefficients iteratively. 
Since our luq and u> can be chosen to be mass- 
independent and the mass dependence of v is very 
weak also on the quantum level for £ > 2 or so 
(see below), we do not expect serious problems. 

The above strategy will require some time and 
effort to complete. Due to the phenomenological 
importance and difficulty of the non-perturbative 
study of heavy quarks we have therefore decided 
to first present a feasibility study of our ap- 
proach that does not tune all coefficients non- 
perturbatively, but only does so for the leading 
one, v, and uses a tree-level tadpole improvement 
estimate || for ujq and u>. For the tadpole esti- 
mate we use the mean links in Landau gauge. For 
the isotropic case it is known [[)) that this gives 
an estimate much closer to the non-perturbative 
value (H,^) than the mean plaquette estimate. 

We performed simulations of the charmonium 
spectrum at several couplings for £ = 3 and £ = 2, 
corresponding to spatial lattice spacings in the 
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o Aniso SW, f=3 
o Aniso SW, f=S 
x Fermilab [ll] 
□ NRQCD v B , u p [12] - 
NRQCD v s , u L [12] 
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Figure 3. Quenched cc hyper-fine splitting l 3 Si — 
l 1 ^ from our (o, o) and other methods. A joint 
a 2 fit of our results is shown. For details see 101. 



range 0.17 to 0.3 fm (the scale was set very accu- 
rately via r JTo|). The bare quark mass mo and 
velocity of light v were tuned so that the spin 
average IS meson mass equals its observed value 
and that c(0) = 1 for the pseudo-scalar. c(0) was 
obtained by extrapolation from c(p) at the two 
lowest on-axis momenta. We find that the v have 
a rather weak mass dependence; they are even 
quite close to their classical values of figure [|. We 
compared the scale from setting r§ = 0.5 fm with 
that from the spin-averaged IP — IS splitting, 
finding that they agree within a few percent. 

Our main results are the hyper-fine splittings 
for the S- and P-states shown in figures || and [|, 
where they are compared with results from other 
approaches. Each of our data points uses 400 — 
600 configurations with suitably smeared quark 
sources for early plateaux. We regard our re- 
sults as the most reliable (note that the strong 
mass dependence of v, ujq and u) is currently not 
taken into account in the heavy relativistic ap- 
proach; also, it had been concluded earlier [Q 
that NRQCD is breaking down for charmonium, 
at least for spin splittings). We emphasize the 
rather non-trivial agreement between the contin- 
uum limits of our results for £ = 2 and 3 in an 
a 2 continuum extrapolation, which suggests that 
our estimates of ujq, lo eliminate most O(a) errors. 

To summarize, on-shell O(a) improvement on 
anisotropic lattices is a computationally effi- 
cient scheme for the study of heavy quarks 
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Figure 4. P-state splitting l 3 Pi - 1 3 P . 

(all our results were obtained on workstations). 
Its efficiency is chiefly due to two factors: (i) 
Anisotropic lattices give very accurate effective 
masses, (ii) Tuning is, which is the correct way to 
proceed anyhow, is in practice simpler and more 
accurate than using the "kinetic mass prescrip- 
tion" as in other approaches. Probably most im- 
portant, however, it seems that systematic errors 
due to mass-dependent, relativistic, and quantum 
corrections to the coefficients in the action can be 
more easily eliminated on anisotropic lattices jl . 
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